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ABSTRACT 


The  yiermal  expansion  of  coupons  of  3D  carbon-carbon  is  analyzed  in  terms  of  the  thermo-elastic 
properties'  of  the  mini-constituents  (yarn  bundles  and  matrix  pockets)  and  the  stress  transfer  across 
their  interfaces.  The  weakness  of  the  interfaces  limits  the  stress  transfer  by  shear.  When  the 
shear  strength  of  the  interface  is  exceeded  and  debonding  occurs,  stress  transfer  in  the  de-bonded 
region  is  assumed  to  be  frictional.  The  frictional  shear  capability  is  estimated  as  a  function  of 
the  inter-constiiuent  stresses  that  occur  on  heating.  The  analysis,  which  deals  with  the 
stress-transfer  using  a  shear-lag  approach,  provides  estimates  of  the  displacements  between 
mini-constituents.  Numerical  results  are  presented  to  show  how  elastic  and  frictional  shear  lag 
influences  measurements  of  thermal  expansion.  The  relation  of  specimen  dimensions  to  the  geometry 
of  the  composite's  unit  cell,  and  details  of  the  measurement  technique,  are  predicted  to  influence 
the  measured  expansions)  for  some  techniques,  the  measured  expansion  is  predicted  to  differ 
significantly  from  the  true  expansion  of  the  bulk  composite  material.  Analytical  predictions  are 
shown  to  compare  reasonably  well  to  Lander's  data  on  end  effects,  up  to  about  A000  F,  beyond  which 
the  likelihood  of  creep/relaxation  phenomena  makes  the  analysis  inapplicable.  ^Directions  for 
improving  the  analysis  are  described,  and  requirements  for  experimental  determnation  of  input 
properties  are  identified.  Recommendations  for  accurate  thermal  expansion  testing  are  provided. 


INTRODUCTION 


Thermal  expansion  is  usually  measured  by  observing  the  change  in  length  of  a  uniformly  heated 
rod  or  bar  of  material.  For  many  materials,  those  having  a  microstructural  dimension  that  is  small 
enough  that  the  material  may  be  viewed  as  homogeneous  on  the  scale  of  the  specimen  rod  or  bar,  the 
length  change  in  the  test  is  a  direct  measure  of  the  thermal  expansion  of  the  bulk  material. 

However,  for  3D  carbon-carbon  composites,  the  cross-sectional  dimension  of  the  yarns  is  in  the  order 
of  1  mm,  which  is  not  very  small  compared  to  the  dimensions  of  typical  specimens.  Insofar  as  the 
composite's  thermal  expansion  depends  on  stress  interactions  among  the  various  yarns  and  the  matrix 
pockets,  and  insofar  as  these  stresses  are  influenced  by  the  presence  of  free  surfaces  at  the 
boundaries  of  the  specimen,  the  specimen's  change  in  length  is  not  necessarily  a  direct  measure  of 
the  thermal  expansion  of  the  bulk  composite. 

The  analysis  presented  in  this  paper  is  intended  to  model  stress  transfer  between  constituents 
of  a  3D  composite  near  free  surfaces  and  its  effects  on  the  length  changes  that  would  be  measured  in 
a  thermal  expansion  test.  The  inputs  to  the  analysis  include  the  thermo-elastic  properties  of  the 
mini-constituents,  as  the  yarn  bundles  and  matrix  pockets  are  sometimes  called,  the  strength  of  the 
interface  between  a  yarn  and  the  rest  of  the  composite,  the  dimensions  of  the  composite's  unit  cell, 
and  the  length  of  the  specimen.  Outputs  include  the  change  in  specimen  length  on  heating,  as  a 
function  of  how  it  is  measured,  and  the  stresses  generated  by  heating  within  the  yarn  and  the 
surrounding  composite.  Recommendations  for  avoiding  large  errors  in  estimating  thermal  expansion  of 
the  composite  can  be  derived  from  the  results. 


SHEAR  LAG  ANALYSIS 


As  shown  in  Figure  1,  the  3D  composite  is  considered  to  consist  of  two  phases!  "yarn", 
comprising  the  primary  bundles  that  are  oriented  axially  (normal  to  the  free  surface),  and  "matrix" 
comprising  the  other  yarn  bundles  and  the  matrix  pockets.  When  the  composite  is  heated,  differences 
in  the  thermal  expansions  of  the  two  phases  give  rise  to  stresses.  Generally,  on  heating  of 
carbon-carbons,  the  yarn  will  be  in  axial  tension  while  the  matrix  will  be  in  axial  compression; 
perpendicular  to  the  axis,  in  the  transverse  direction,  the  yarn  will  be  in  compression  (Ref.  1  and 
2,  for  example).  At  the  free  surface,  in  the  absence  of  externally  applied  tractions,  the  stresses 
in  yarn  and  matrix  will  be  zero.  Thus,  there  is  a  region  near  the  surface  in  which  the  stresses 
vary,  which  implies  the  existence  of  shear  at  the  interface  between  yarn  and  matrix.  The  analysis 
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Fig.  1.  Idealization  of  3D  composite. 


Fig.  2.  Unit  cell  of  specimen. 


described  below  defines  the  stress  gradients  in  terms  of  a  simple  shear  lag  model  (adapted  from  Ref. 
3).  Once  the  stresses  are  defined,  the  strains  and  displacements  can  be  calculated  for  yarn  and 
matrix. 


Specifically,  we  consider  simple  rectangular  bars  such  as  are  commonly  used  for  thermal 
expansion  measurements.  The  symmetry  of  the  situation  allows  us  to  do  the  analysis  by  studying  half 
the  length  of  such  a  specimen,  considering  a  single  yarn  and  its  surrounding  matrix  (Fig.  2). 

Because  of  the  shear  lag,  the  matrix  will  displace  axially  more  than  the  yarn,  giving  rise  to  a  wavy 
surface  at  the  specimen  end.  This  "roughening"  of  the  surface  can  produce  discrepancies  between  the 
measured  growth  of  a  specimen  and  the  true  thermal  expansion  of  the  composite  (eg,  Ref.  4). 

If  the  interfacial  shear  stress  exceeds  the  strength  of  the  interface,  debonding  will  occur. 

In  the  debonded  region,  a  frictional  shear  stress  can  exist,  which  depends  in  magnitude  on  the 
compressive  stress  across  the  interface,  among  other  factors.  In  the  region  that  remains  bonded, 
the  shear  stress  is  predicted  from  the  elastic  shear  lag  analysis.  Figures  3  shows  schematically 
the  expected  relative  displacements  of  yarn  and  matrix,  in  the  bonded  and  debonded  regions. 

We  assume  the  yarns  have  square  cross-sections i  it  would  not  be  difficult  to  extend  the 
analysis  to  rectangular-section  yarn  bundles,  but  the  major  points  to  be  made  in  this  paper  can  be 
illustrated  with  the  simpler  analysis  of  square-section  reinforcements.  Many  of  the  equations  for 
circular-section  reinforcements  are  identical  to  those  for  square-section  yarns;  thus,  much  of  the 
analysis  is  also  applicable  to  shear  lag  between  circular  fibers  and  matrix  within  a  yarn  bundle. 


Shear  Lag  Equations 


For  a  specimen  that  is  uniformly  heated  without  any  external  forces,  stress  equilibrium 
requires  the  axial  forces  in  yarn  and  matrix  to  sum  to  zero  at  every  axial  location! 

O'f '  0  (0 

We  assume,  in  accordance  with  the  level  of  this  analysis,  that  and  C„do  not  vary  with  radial 
position. 


Consideration  of  equilibrium,  between  shear  stress  and  axial  stress  in  the  fiber,  givesi 

I  r 


-ir  ri 


rsr. 


CO 


262 


a)  fully  clastic 


b)  brsoNtED  with  friction  : 
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Fig.  3.  Schematic  of  deformations  at  the  free  surface. 


For  the  matrix,  equilibrium  requires: 
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When  the  interface  is  intact,  and  the  yarn  and  matrix  behave  elastically,  we  may  define  an  average 
shear  strain:  ° 
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which  can  be  related  by  Hooke's  law,  to  the  interface  shear  stress: 
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Here,  C  is  an  "effective"  shear  modulus  that  represents  the  fiber  and  matrix  in  appropriate 
proportions  and  accounts  for  our  use  of  maximum  shear  stress  in  combination  with  average  shear 
strain.  This  effective  shear  modulus  may  be  derived  by  expressing  the  average  shear  strain  at  given 
axial  location  in  terms  of  the  average  shear  strains  in  yarn  and  matrix: 
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Applying  Hooke's  law: 
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where  the  average  shear  stresses  are: 
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Introducing  Eq  2  and  A,  and  integrating  Eq  10  and  11,  gives i 
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The  needed  expression  for  effective  shear  modulus  may  now  be  obtained  by  substituting  Equations  12, 
13,  9,  and  8  into  Eq  7i 
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Usually,  G  is  about  1.5  to  3  times  the  low  bound  (uniform  microstress)  estimate  of  the  composite 
shear  modulus,  and  may  in  some  cases  be  larger  than  the  upper  bound  (uniform  microstrain)  estimate. 

To  obtain  a  solution  for  interface  shear  stress  as  a  function  of  axial  position,  x,  we  proceed 
by  differentiating  Eq.  7s 
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Substituting  Eq 

16  and  1  into 

Eq  15, 

we  get: 
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The  maximum  yarn  stress  that  can  be  generated  by  heating  alone  will  occur  at  the  midlength  of  a  very 
long  composite  body,  where  “r^K  will  be  essentially  zero  at  midlength.  From  Eq  ‘7,  the  maximum  yarn 
stress  s*  would  be: 
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Substituting  Eq  19  into  Eq  17,  and  differentiating  again,  we  get: 
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Substituting  the  equilibrium  condition  of  Eq  5  gives  us  a  differential  equation  in 
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A  general  solution  to  Eq  21  is: 
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where  A  and  B  are  constants  that  depend  on  the  boundary  conditions.  We  can  readily  show  B  =  0 
because  the  interface  shear  stress  TJ  must  be  zero  at  x  :  0  (because  of  symmetry  of  the  .«pecimen). 
Thus: 
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For  a  specimen  of  finite  length,  the  gradient^rj4< at  midlength  will  not,  in  general,  be 
the  actual  maximum  yarn  stress  (at  x  =  0)  can  be  expressed  as: 
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Equating  Eq  25  and  17,  at  x  =  0,  and  using  Eq  26,  we  get  a  definition  of  A: 

At  this  point,  is  undefined.  To  proceed,  we  return  to  consideration  of  the  axial  stress  in  the 
yarn.  From  Eq  5  and  26: 
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Introducing  Eq  24  and  integrating: 
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The  yarn  axial  stress  <Sjp  at  the  end  of  the  elastic  region  (x  s  x • )  may  be  estimated  and  used  in  Eq 
29  to  provide  a  second  relation  between  A  and 

A  r  (*) 

If  no  debonding  occurs  in  the  specimen,  and  it  behaves  elastically  over  its  entire  length,  the  fiber 
stress  at  the  end  of  the  elastic  region  (x  =  L)  is  zero.  If  debonding  does  occur,  the  fiber  stress 
at  the  end  of  the  bonded  region  (x  =  x'<  L)  is  given  by  Eq  41  (derived  in  the  next  Subsection). 


We  solve  for  ^  by  equating  Eq  27  and  Eq  30: 
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In  this  "analysis,  we  solve  for  x'  numerically  by  requiring  that  the  shear  stress  at  x'  be  the 
interfacial  shear  strength.  An  iterative  procedure  first  assumes  an  arbitrary  value  of  x',  solves 
for  and  A,  computes  the  elastic  shear  stress  at  x'  from  Eq  24,  and  if  that  value  does  not  equal 
the  shear  strength  (from  Eq  49,  below),  assumes  a  new  value  of  x'  and  repeats  the  calculation  and  so 
on  until  the  discrepancy  is  negligible  (less  than  one  percent  of  the  shear  strength). 


The  change  in  length  of  the  matrix,  in  the  bonded  region,  is: 
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Introducing  Eq  1  and  29,  and  performing  the  integration,  we  get: 
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The  corresponding  length  change  of  yarn,  in  the  bonded  region,  may,  with  the  aid  of  Eq  6,  be 
estimated  as:  ,  . 
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Frictional  Shear  Stress 


If  the  elastic  shear  stress  at  the  primary  yarn  interface  exceeds  the  interface  shear  strength, 
debonding  will  occur  and  the  yarn  may  slip  relative  to  the  matrix.  A  frictional  shear  stress  will 
exist  at  the  interface.  We  assume  this  frictional  shear  stress  T 1  is  the  sum  of  two  factors,  the 
product  of  a  friction  coefficient^M.  and  the  compressive  stress  mrc  acting  across  the  interface, 
and  a  constant  T*  representing  resistance  by  other  effects,  such  as  mechanical  interlocking  between 
rough  interface  surfaces: 


The  transverse  compressive  stress  arises  from  the  minimechanical  interactions  between  the  transverse 
yarns  and  the  rest  of  the  composite  (which  includes  our  primary  axial  yarn).  Because  shear  lag 
phenomena  apply  also  to  the  transverse  yarns,  the  compression  will  vary  with  distance  from  a 
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transverse  free  surface.  At  a  transverse  surface,  the  compressive  stress  on  the  primary  yarn  will 
approach  zero.  Toward  the  center  of  a  large  enough  body,  the  compressive  stress  will  approach  a 
maximum  value  that,  for  each  set  of  transverse  yarns  may  be  estimated  as  (eg,  Ref  2 ) t 
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where  3%  is  the  composite's  thermal  expansion: 
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and  vj  is  the  yarn  volume  fraction  in  the  appropriate  direction.  For  illustrative  purposes,  we 
simplify  the  analysis  to  deal  only  with  square  primary  yarns  in  a  matrix  characterized  by  equal 
distribution  of  yarns  in  the  two  transverse  directions.  Then: 
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and  vT  is  the  total  volume  fraction  of  yarn  in  the  composite,  usually  in  the  vicinity  of  .7  to  .75. 
The  task  of  estimating  the  compressive  stress  distribution  in  a  slender  bar  specimen  is  not  trivial 
as  it  would  require  iterative  application  of  the  shear  lag  analysis  to  the  primary  and  transverse 
yarns.  To  illustrate  the  effects,  we  define  a  factor  : 
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and  we  consider  two  examples  in  the  application  of  the  analysis:  1)  <4>  =  0.1  as  applying 
approximately  to  the  corners  of  a  specimen,  and  2)  If'  =1.0  (maximum  compression)  as  providing  a 
bound  to  the  situation  at  the  specimen  centerline  (Figure  A). 

We  assume  further  that  the  frictional  shear  stress  is  invariant  with  axial  position.  Then,  by 
integrating  Eq.  5,  we  find  the  distribution  of  axial  stress  in  the  yarn  to  be: 


(«') 


The  change  in  length  of  the  yarn  in  the  debonded  region  is: 
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and  the  length  change  of  matrix  is: 
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Performing  the  integrations,  we  get: 
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and 
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Fig.  4.  Yarn  positions. 


Total  Length  Change  of  the  Specimen  on  Heating 


From  Eq  44  and  45,  and  Eq  33  and  34,  it  is  clear  that  the  length  change  measured  in  a  thermal 
expansion  test  will  depend  on  whether  the  test  technique  reads  the  length  change  of  the  matrix  or  of 
the  yarn.  The  total  length  change  of  the  matrix  region  is  the  sum  of  the  length  changes  in  the 
bonded  and  debonded  regions: 
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Similarly,  the  total  length  change  of  the  yarn  is: 
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The  difference  between  these  two  length  changes  is  the  surface  roughening  shown  in  Figure  3. 
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INPUT  PROPERTIES 


The  illustrative  results  presented  later  are  intended  to  represent  graphitized  3D  carbon-carbon 
material  made  with  T-300  yarns  with  pitch  densification,  such  as  used  in  the  7-inch  and  15-inch 
billet  programs  (Ref.  5  and  6).  This  section  describes  the  approach  to  estimating  the  material 
properties  used  in  the  analysis.  The  values  used  to  produce  the  illustrative  results  (next  Section) 
are  listed  in  Table  1. 

The  axial  modulus  and  thermal  expansion  of  the  yarn  phase  are  derived  from  test  data.  For 
modulus,  the  correlation  by  F.  I.  Clayton,  reported  in  Ref.  6,  provides  estimates  of  in-situ 
filament  modulus  Efil  vs.  temperature.  We  assume  the  yarn  is  composed  of  60  percent  filament  by 
volume,  so  that  Ef  5  .6  x  Efil.  For  axial  thermal  expansion  of  the  yarn,  we  use  directly  the  data 
obtained  by  Lander  (Ref.  7)  on  a  specimen  comprising  a  single  yarn  bundle. 

The  transverse  expansion  of  the  matrix  phase  (that  includes  the  transverse  yarns  and  matrix 
pockets)  is  based  directly  on  the  measurements  by  Lander  (Ref.  7)  on  a  specimen  from  which  the  axial 
yarns  were  excised. 

The  remaining  properties  of  the  matrix  phase  and  of  the  yarn  are  estimated  using  simple 
rule-of-mixtures  equations  from  values  for  filaments  and  pitch-precursor  matrix  recommended  by 
Kibler  (eg,  Ref.  8).  Estimates  are  required  here  because  direct  experimental  data  are  unavailable. 

Graphitized  3D  carbon-carbons  are  known  to  have  a  regularly  cracked  microstructure  iRef.  2,  for 
example).  In  the  current  context,  the  existence  of  these  microcracks  affects  the  thermo-elastic 
properties  of  the  matrix  region,  the  transverse  elastic  properties  of  the  yarn,  and  the  effective 
shear  modulus  used  in  the  shear  lag  analysis.  The  simplest  approach  to  accounting  for  the 
microcracks  is  to  assign  efficiency  factors  to  the  extensional  and  shear  moduli  of  the  matrix 
material  (eg,  Ref.  8  and  9)s 


where  the  "starred"  properties  represent  uncracked  material.  Following  the  recommendation  of  Ref. 

9,  we  assign  a  lower  efficiency  to  the  extensional  behavior  than  to  the  shear  behavior.  From 
Kibler's  work,  we  use  a  shear  efficiency  factor  of  O.A.  From  attempts  to  match  some  of  Lander's 
thermal  expansion  data  (Ref.  7),  it  appears  that  the  extensional  efficiency  factor  is  in  the 
vicinity  of  .05  to  .10. 

Data  for  interfacial  shear  strength  and  friction  coefficients  is  lacking.  From  a  brief  review 
of  related  information  (eg,  pull-out  failures  in  tensile  tests,  microshear  data,  shear  tests  of  ID 
and  2D  composites,  and  friction  coefficients  measured  at  room  temperature  between  composite  parts) 
approximate  values  were  estimated.  The  values  used  in  the  illustrative  calculations  are  listed  in 
Table  1. 

The  input  value  for  shear  strength,  1gL  ,  represents  the  strength  in  the  absence  of  substantial 
compression  across  the  interface.  For  high  values  of  flfe  it  is  possible  that  the  estimated  friction 
shear  stress  (Eq  35)  will  be  greater  than  the  input  shear  strength.  While  other  approaches  are  also 
attractive  (eg,  use  of  a  multiaxial  fracture  criterion  for  shear  strength),  in  this  paper  we  take 
the  effective  elastic  shear  strength,  ,  to  be  the  larger  of  the  input  strength  and  the  estimated 
friction  shear  stress: 

rei  *  ***  C1**#  r'3  G«) 


TABLE  1.  INPUTS  TO  ILLUSTRATIVE  ANALYSIS 

VALUE  AT  TEMPERATURE 


PROPERTY 

Temp,  F 

0 

1000 

2000 

3000 

A000 

shear  modulus,  yarn 

psi 

Gf* 

1.50E+06 

1.51E+06 

1.60E+06 

1.71E+06 

1.A9E+06 

shear  mod,  transverse  composite 

psi 

Gm# 

5.36E+05 

5.OE+05 

5.75E+05 

6.12E+05 

5.3AE+05 

Young's  modulus,  yarn,  axial 

psi 

Ef 

3.60E+07 

3.60E+07 

3.60E+07 

2.A0E+07 

2.10E+07 

Young's  modulus,  yarn,  transverse 

psi 

EfT" 

1 . 38E+06 

1.A1E+06 

1.50E+06 

1.61E+06 

1.A2E+06 

Young's  modulus,  transv.  composite 

psi 

Em** 

1 . 38E+06 

1.A1E+06 

1.50E+06 

1.61E+06 

1.A2E+06 

extensional  bond  efficiency 

etaE 

0.05 

0.05 

0.05 

0.05 

0.05 

shear  bond  efficiency 

etaS 

0.A0 

0.  A0 

0.  A0 

0.  A0 

0.  A0 

CTE,  yarn,  axial 

per  F 

alpf 

0.00E+00 

-2.00E-07 

2.00E-07 

A.00E-07 

5.50E-07 

CTE,  yarn,  transverse 

per  F 

alpfT 

0. OOE+OO 

A.6AE-06 

A.A6E-06 

A.33E-06 

A.89E-06 

CTE,  transverse  composite 

per  F 

a  1pm 

0. 00E+00 

1.60E-06 

2.25E-06 

3.00E-06 

3.75E-06 

nominal  interface  shear  strength 

psi 

tauEL* 

1000 

1000 

1000 

1000 

1000 

minimum  frictional  shear  strength 

psi 

tau' 

200 

200 

200 

200 

200 

coeff.  friction 

mu 

0.3 

0.3 

0.3 

0.3 

0.3 

total  yarn  volume  fraction 

vY 

0.75 

0.750 

0.750 

0.750 

0.750 

ILLUSTRATIVE  RESULTS 


The  shear  lag  analysis  Is  applied  here  to  the  prediction  of  thermal  expansion  data  as  a 
function  of  the  technique  used  for  measurement  of  specimen  length  change.  The  various  specimen 
configurations  and  length  measurements  used  for  3D  carbon-carbons  are  sketched  in  Figure  5. 

The  first  example  deals  with  analysis  of  "matrix"  expansion  u„f  and  yarn  expansion  u^#, 
providing  calculations  that  may  be  compared  to  Lander's  data  (Ref.  7)  for  flat-ended  and 
protruding-yarn  specimens,  respectively.  Lander  tested  2-inch  long  specimens  taken  from  the 
outer-diameter  region  of  a  woven  cylindrical  billet,  oriented  in  the  billet's  axial  direction.  We 
used  the  input  properties  shown  in  Table  1  and  implemented  the  analysis  using  an  estimated  yarn 
volume  fraction  of  0.19  and  a  unit-cell  radial  dimension,  r„,  of  .07  inch.  Figure  6  shows  the 
analytical  predictions  to  be  in  reasonable  agreement  with  Lander's  data  up  to  about  2200  C  (4000  F). 
Correlations  above  4000  F  were  not  attempted  because  the  analysis  does  not  currently  include 
creep/relaxation  effects,  which  are  important  at  higher  temperatures. 

The  predictions  are  shown  as  lines  in  Figure  6|  the  top  solid  line  is  the  prediction  for  the 
corner  of  a  flat  specimen  where  the  frictional  stress  is  reduced  because  minimal  compression  can 
exist  across  the  yarn  interface;  the  bottom  solid  line  is  the  prediction  for  the  centerline  region 
of  a  flat  specimen  that  is  sufficiently  large  in  cross-section  to  exhibit  maximum  compression  at  the 
yarn  interface.  The  lines  tend  to  bracket  Lander's  data,  suggesting  that  the  analytical  model  is 
capable  of  predicting  the  major  effects  in  such  specimens.  Also  shown  is  the  predicted  response  of 
the  yarn,  which  corresponds  to  data  from  protruding-yarn  specimens.  The  difference  between 
predicted  responses  of  yarns  at  the  corner  and  centerline  is  negligible,  so  only  one  dotted  line  is 
plotted.  In  Figure  7,  these  predictions  are  compared  to  the  theoretical  response  of  the  bulk 
composite  (estimated  by  Eq  37);  we  see  that  the  yarn  responses  are  quite  close  to  the  theoretical 
value,  whereas  the  matrix  response  is  a  substantial  overestimate. 
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Pig.  5.  Various  end  configurations  used  in  thermal  expansion  testing. 


Fig.  6.  Comparison  of  analysis  tc  Lander's  data  for  flat  and  protruding-yarn  specimens. 
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Fig.  7.  Discrepancies  between  predicted  measurements  and  theoretical  expansion  of  bulk  composite. 


Figure  8  shows  the  predicted  roughening  (in  hundredths  of  an  inch)  of  a  flat  end-surface  and 
the  predicted  extent  of  debonding.  Figure  9  shows  the  predicted  variation  of  yarn  axial  stress  and 
interface  shear  stress,  for  the  centerline  region. 

Additional  analyses  were  done  to  explore  the  effects  of  specimen  length  (actually,  the  key 
variable  is  the  ratio  of  specimen  length  to  unit-cell  dimension,  r„)  and  yarn  volume  fraction.  For 
these  calculations,  we  assumed  the  yarn  experiences  the  maximum  transverse  compression.  Thus  the 
data  provide  a  low-bound  expansion  for  the  centerline  of  flat-ended  or  spherical-ended  specimens. 
Also,  the  interface  shear  strength  was  taken  as  2000  psi,  rather  than  the  1000  psi  used  in  analyzing 
Lander's  data. 

The  effect  of  specimen  length  on  predicted  expansion  data  to  4000  F  is  shown  in  Figure  10,  for 
a  yarn  fraction  of  .21  and  a  unit-cell  radius  of  .07-inch.  We  see  that  the  potential  errors  from 
using  flat-ended  specimens  decrease  as  specimen  length  increases.  It  may  be  of  interest  to  note 
that  the  specimen  lengths  used  for  expansion  measurements  on  fine-weave  carbon-carbons  (eg,  Ref.  10) 
were  sufficient  to  make  such  errors  nearly  negligible;  to  achieve  similar  accuracy  with  the  coarse 
weave  composites  such  as  the  billet  Lander  studied,  we  would  need  specimen  lengths  in  the  order  of 
ten  inches.  The  predicted  surface  roughening  and  the  length  of  debond  are  insensitive  to  specimen 
length,  in  the  range  studied.  The  yarn  stress  generated  at  the  specimen  midlength  is  within  one 
percent  of  the  theoretical  value,  even  for  the  shortest  length  shown;  this  is  consistent  with  the 
length  of  the  shear-lag  region  shown  in  Fig  9.  Thus,  the  assumption  of  full  compression 

across  the  centerline  yarn  interface  may  be  close  to  reality  for  specimens  about  3/8-inch  in 
cross-section.  However,  as  noted  earlier,  the  proper  analysis  of  the  transverse  compression  is 
complex,  so  the  implication  Just  stated  should  viewed  cautiously. 

The  effects  of  yarn  volume  fraction  on  predicted  expansion  data  to  4000  F  are  shown  in  Figure 
11,  for  2-inch  long  specimens  having  a  unit-cell  radius  of  .07  inch.  Potential  errors  in  the  use  of 
flat-ended  specimens  increase  significantly  at  low  values  of  yarn  volume  fraction.  This  finding  may 
aid  in  explaining  some  of  the  very  high  values  of  radial-direction  expansion  reported  for  some 
cylindrical  billet  materials  (eg,  Ref.  6). 
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DISTANCE  FROM  SPECIMEN  CENTER 

Fig.  9.  Predicted  stress  gradients  at  interface  and  in  yarn  at  4000  F  at  specimen  centerline 


PERCENT  ERROR 


Fig.  10.  Effect  of  specimen  half-length  on  predicted  errors  at  4000  F. 


Fig.  11.  Effect  of  yarn  volume  fraction  on  predicted  errors  at  4000  F. 


CONCLUDING  REMARKS 


The  shear  lag  analysis  shows  that  thermal  expansion  data  can  be  sensitive  to  the  details  of  the 
measurement  technique,  especially  to  the  length  of  specimen  used  and  to  the  manner  of  measuring  its 
change  in  length  (Figure  5).  Potential  errors  are  larger  for  specimens  of  smaller  volume  fraction 
of  yarns  in  the  direction  of  measurement,  and  for  coarse-weave  composites. 

On  the  basis  of  the  findings,  the  use  of  flat-ended  specimens  should  be  discouraged.  The  use 
of  spherical-end  specimens  is  preferable  to  flat-end  specimens)  however,  significant  error  may  be 
experienced  due  to  surface  roughening.  The  use  of  pin-ended  specimens  should  provide  quite  accurate 
data  if  the  pin  rests  on  the  end  of  a  yarn  and  is  of  a  diameter  smaller  than  the  yarn  cross-section; 
otherwise,  data  from  pinned  specimens  may  be  influenced  by  roughening  at  the  base  of  the  pin.  (The 
shear  lag  analysis  could  readily  be  extended  to  treat  the  effect  of  pin  depth).  Of  all  the 
configurations  shown  in  Figure  5,  the  most  accurate  appears  to  be  the  protruding-yarn  specimen 
(Fig.5c/4  which  was  first  used  by  Lander  in  the  end-effects  study  of  Ref.  7. 

Many  of  the  properties  that  are  inputs  to  the  analysis  are  not  well  known.  It  is  recommended 
that  research  be  directed  toward  measuring  transverse  properties  of  yarn  bundles,  in-situ  properties 
of  matrix-pocket  regions,  establishing  the  appropriate  efficiency  factors  to  account  for 
pre-existing  microcracks,  measuring  yarn  interface  strengths  and  friction  coefficients,  all  as 
functions  of  temperature.  While  some  data  of  the  types  listed  is  available,  more  is  needed  for 
specific  materials  of  interest. 

It  would  be  worthwhile  to  attempt  experimental  verification  of  the  predicted  debond  lengths  and 
surface  roughenings. 

Currently,  the  analysis  does  not  treat  creep/relaxation  effects,  which  are  undoubtably 
important  at  temperatures  above  4000  F.  Extension  of  the  analysis,  and  the  acquisition  of  relevant 
constituent  creep  data  would  be  worthwhile.  Other  extensions  of  the  analysis  should  also  be 
attempted.  These  include  treatment  of  rectangular-section  yarns,  allowance  for  differing  volume 
fractions  for  the  two  sets  of  transverse  yarns,  and  better  treatment  of  the  compressive  stresses 
that  influence  the  frictional  shear. 

Motivation  for  the  suggested  further  work  includes  the  belief  that  shear  lag  analysis  is  a 
powerful  tool  for  studying  the  structural  response  of  3D  carbon-carbon  components  near  free 
surfaces.  In  comparison  to  finite-element  analyses  of  mini-mechanical  interactions,  the  shear  lag 
approach  is  simple  and  inexpensive)  the  simplifications  (including  neglect  of  Poisson's 
interactions)  and  lack  of  rigor  (in  the  sense  of  elasticity  solutions)  are  compensated  for  by  the 
ability  to  treat  approximately  such  phenomena  as  inter-constituent  slip,  which  are  difficult  to 
model  .-igorously. 


NOMENCLATURE 

oi  =  thermal  expansion  coefficient,  secant  value 
%  =  shear  strain 
9  s  extensional  stress 
T  =  shea.’  stress 
E  =  Young's  modulus 
G  =  shear  modulus 

L  =  half  the  length  of  the  specimen 
S  s  ratio  of  yarn  radius  to  matrix  radius 
T  -  temperature 
e  =  extensional  strain,  total 
r  =  radial  distance,  or  radius 
u  =  axial  displacement 
v  =  volume  fraction  of  yarn 
x  s  ax'il  distance 
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Subscripts 

c  =  compression 

f  =  pertaining  to  yarn 

i  =  pertaining  to  interface 

j  =  index  defining  a  yarn  direction 

m  =  pertaining  to  matrix 

T  =  transverse  to  the  yarn  axis 
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